A high-sensitivity laser-pumped Mj-magnetometer 
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We discuss the design and performance of a laser-pumped cesium vapor magnetometer in the M^- 
configuration. The device will be employed in the control and stabilization of fluctuating magnetic 
fields and gradients in a new experiment searching for a permanent electric dipole moment of the 
neutron. We have determined the intrinsic sensitivity of the device to be 29 IT in a 1 Hz bandwidth, 
limited by low frequency noise of the laser power. In the shot noise limit the magnetometer can 
reach a sensitivity of 10 IT in a 1 Hz bandwidth. We have used the device to study the fluctuations 
of a stable magnetic field in a multi-layer magnetic shield for integration times in the range of 2-100 
seconds. The residual fiuctuations for times up to a few minutes are traced back to the instability 
of the power supply used to generate the field. 

PACS numbers: 07.55.Ge, 32.30.Dx, 32.70.Jz, 33.55.-b 
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I. INTRODUCTION 

In many areas of fundamental and applied science the 
sensitive detection of weak magnetic fields and small field 
fluctuations is of great importance. In the applied sec- 
tor this concerns, for instance, non-destructive testing of 
materials 1], geomagnetic and archaeological prospect- 
ing j2,] , and the expanding field of biomagnetism Q ■ I^i 
the realm of fundamental physics, strong demands on 
magnetometric sensitivity are placed by modern experi- 
ments looking for small violations of discrete symmetries 
in atoms and elementary particles. For instance, many 
experiments searching for time-reversal or parity viola- 
tion rely on the precise monitoring and control of mag- 
netic fields, with the sensitivity of the overall experiment 
directly related to the ultimate sensitivity and stability of 
the magnetic field detection. Picotesla or even femtotesla 
sensitivity requirements for averaging times of seconds to 
minutes are common in that field. 

Our particular interest in this respect lies in the search 
for a permanent electric dipole moment (EDM) of the 
neutron. Such a moment violates both time reversal in- 
variance and parity conservation. A finite sized EDM 
would seriously restrict theoretical models that extend 
beyond the standard model of particle physics 0| . Re- 
cently our team has joined a collaboration aiming at a 
new measurement of the permanent EDM of ultra-cold 
neutrons (UCN) to be produced from the UCN source 
under construction at Paul-Scherrer-Institute in Switzer- 
land [3- A neutron EDM spectrometer will be used, 
in which the neutron spin-flip frequency will be mea- 
sured by a Ramsey resonance method in UCN storage 
chambers exposed to a homogenous magnetic field. Each 
neutron chamber has two compartments in which the 
neutrons are exposed to a static electric field oriented 
parallel/antiparallel to the magnetic field. The signa- 
ture of a finite EDM will be a change of the neutron 
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Larmor frequency that is synchronous with the reversal 
of the relative orientations of the magnetic and electric 
fields. Magnetic field instabilities and inhomogeneities 
may mimic the existence of a finite neutron EDM. The 
control of such systematic effects is therefore a crucial fea- 
ture of the EDM experiment. It is planned to use a set of 
optically pumped cesium vapor magnetometers (0PM), 
operated in the M^ configuration [a, l3 ^^ perform that 
control. 

Although OPMs pumped by spectral discharge lamps 
are suited for the task, we have opted for a system of laser 
pumped OPMs (LsOPM). It was shown previously that 
the replacement of the lamp in an 0PM by a resonant 
laser can lead to an appreciable gain in magnetometric 
sensitivity j3, l3l- Laser pumping further offers the ad- 
vantage that a single light source can be used for the 
simultaneous operation of several dozens of magnetome- 
ter heads. In that spirit we have designed and tested 
a LsOPM with a geometry compatible with the neu- 
tron EDM experiment. In this report we present the 
design and the performance of the Cs-LsOPM operated 
in a phase-stabilized mode and discuss a systematic effect 
specifically related to laser pumping. 



II. THE OPTICALLY-PUMPED M^ 
MAGNETOMETER 

Optically pumped magnetometers can reach extreme 
sensitivities of a few fl/^/ilz |3, comparable to stan- 
dard SQUID (superconducting quantum interference de- 
vice) detectors. Recently a low field 0PM with a sub-fT 
resolution was demonstrated Q- The use of OPMs for 
the detection of biomagnetic signals was recently demon- 
strated by our group |10l[l]|. 

As a general rule the optimum choice of the 0PM 
depends on the specific demands (sensitivity, accuracy, 
stability, bandwidth, spatial resolution, dynamic range, 
etc.) of the magnetometric problem under considera- 
tion. In our particular case the main requirements are 
a highest possible sensitivity and stability for averaging 



times ranging from seconds up to 1000 seconds in a 2 /iT 
field together with geometrical constraints imposed by 
the neutron EDM experiment. 

Optically pumped alkali vapor magnetometers rely on 
an optical-radio frequency (r.f.) resonance technique and 
are described, e.g., in p. When an alkali vapor is irradi- 
ated with circularly polarized light resonant with the Di 
absorption line (transition from the nSi/2 ground state 
to the first nPi/2 excited state), the sample is optically 
pumped and becomes spin polarized (magnetized) along 
the direction of the pumping light. While lamp pumped 
OPMs simultaneously pump all hyperfine transitions of 
the Di line, the use of a monomode laser in a LsOPM 
allows one to resolve the individual hyperfine transitions 
provided that their Doppler width does not exceed the 
hyperfine splitting in both the excited and the ground 
states. This is, for example, the case for the Di tran- 
sition of the alkali isotopes ^^'^Cs and *^Rb. In that 
case it is advantageous to set the laser frequency to the 
F = I + 1/2 -^ F = I — 1/2 transition, which allows one 
to optically pump the atoms into the two (non-absorbing) 
dark states |nS'i/2; F', Mp = F, F — 1) using ct+ polarized 
radiation. A magnetic field Bi(i) oscillating at the fre- 
quency Wrf , which is resonant with the Zeeman splitting 
of the states, drives population out of the dark states into 
absorbing states, so that the magnetic resonance transi- 
tion can be detected via a change of the optical transmis- 
sion of the vapor. That is the very essence of optically 
detected magnetic resonance. 

In the so-called Mx or 45° configuration the static mag- 
netic field Bq to be measured is oriented at 45° with re- 
spect to the laser beam, while the oscillating magnetic 
field Bi(t) is at right angles with respect to Bq (Fig. ^). 
In classical terms, the Larmor precession of the magne- 
tization around Bq (at the frequency lul) is driven by 
the co-rotating component of the Bi(i)-field, which im- 
poses a phase on the precessing spins. The projection of 
the precessing polarization onto the propagation direc- 
tion of the light beam then leads to an oscillating mag- 
netization component along that axis, and therefore to 
a periodic modulation of the optical absorption coeffi- 
cient. The system behaves like a classical oscillator, in 
which the amplitude and the phase of the response (cur- 
rent from a photodiode detecting the transmitted laser 
intensity) depend in a resonant way on the frequency of 
the Bi field. From the resonance condition uj^ — Wrf the 
Larmor frequency and hence the magnetic field can be 
inferred. 

When the AC component of the detected optical signal 
is transmitted to the coils producing the Bi(^) field with 
a 90° phase shift and an appropriate gain, the system will 
spontaneously oscillate at the resonance frequency. In 
that self-oscillating configuration the 0PM can in princi- 
ple follow changes of the magnetic field instantaneously 
with a bandwidth limited by the Larmor frequency only 

Here we have used an alternative mode of operation, 
the so-called phase-stabilized mode. The in-phase ampli- 



tude X, the quadrature amplitude Y and the phase (j) of 
the photocurrent with respect to the oscillating magnetic 
field are given by 
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where x — [vl — t'rf)/Ai4jvvHM is the detuning normal- 
ized to the (light-power dependent) half width at half 
maximum Ai^hwhm of the resonance. 5' is a saturation 
parameter which describes the r.f. power broadening of 
the line. It is interesting to note that the width of the 
phase dependence, which is determined by the ratio of 
the X{x) and Y{x) signals, is independent of S, and 
hence immune to r.f. power broadening. The phase (j){x) 
changes from 0° to —180° as v^t is tuned over the Lar- 
mor frequency. Near resonance the phase is —90° and 
has a linear dependence on the detuning v^ — ^'rf ■ </'(a;) is 
detected by a phase sensitive amplifier (lock-in detector) 
whose phase output drives a voltage-controlled oscillator 
(VCO) which feeds the r.f. coils. The VCO signal, phase 
shifted by 90°, serves as a reference to the phase detector. 
This feedback loop thus actively locks the r.f. frequency 
to the Larmor frequency and the magnetometer tracks 
magnetic field changes in a phase coherent manner. That 
mode of operation is a modification of the self-oscillating 
magnetometer in the sense that the lock-in amplifier, the 
loop filter (PID), and the VCO represent the components 
of a tracking filter which shifts the detected signal by 90° 
and apphes the filtered signal to the r.f. coils. The dif- 
ferences to the self-oscillating scheme are the following: 
the bandwidth of the phase-stabilized magnetometer is 
determined by the transmission function of the feedback 
loop, and the phase shift is always 90° independent of the 
Larmor frequency, while in the self-oscillating scheme the 
phase-shifter has a frequency dependence and is 90° only 
for a single Larmor frequency. Note that the tracking 
filter in a strict sense is not a phase-locked loop (PLL), 
since there is only one detectable frequency in the sys- 
tem, i.e., v-cf. A detuning between the r.f. frequency 
and the Larmor frequency produces a static phase shift, 
while in a real PLL the detuning between the reference 
frequency and the frequency which is locked produces a 
time dependent phase shift. 



III. MAGNETOMETER HARDWARE 

The LsOPM for the n-EDM experiment consists of two 
parts: a sensor head containing no metallic parts except 
the r.f. coils, and a base station mounted in a portable 
19" rack drawer, which contains the frequency stabilized 
laser and the photodetector. The laser light is carried 
from the base station to the sensor head by a 10 m long 
multimode fiber with a core diameter of 800 /im. The 




FIG. 1: Schematic setup of the phase-stabilized magnetome- 
ter in the closed-loop (A) and the scanning (B) mode. The 
dashed box indicates the sensor head. L: lens, P: polariz- 
ing beamsplitter, A/4: quarter-wave plate, PD: photodiode, 
VRG: voltage ramp generator, VCO: voltage-controlled oscil- 
lator, PID: feedback amplifier. The stabilization system of 
the laser frequency is not shown. 



light transmitted through the cell is carried back to the 
detection unit by a similar fiber. The sensor head is de- 
signed to fit into a tube of 104 mm diameter, coaxial with 
the 2 /xT field, and has a total length of 242 mm. The 
main component of the sensor is an evacuated glass cell 
with a diameter of 7 cm containing a droplet of cesium 
in a sidearm connected to the main volume. A constric- 
tion in the sidearm minimizes the collision rate of vapor 
atoms with the cesium metal. The probability of spin 
depolarization due to wall collisions with the inner sur- 
face of the glass cell is strongly reduced by a thin layer 
of paraffin coating the cell walls. The cell was purchased 
from MAGTECH Ltd., St. Petersburg, Russia. A pair 
of circular coils (70 mm diameter separated by 52 mm) 
encloses the cell and produces the oscillating magnetic 
field Bi(i). 

The light driving the magnetometer is produced by 
a tunable extended-cavity diode laser in Littman con- 
figuration (Sacher Lasertechnik GmbH, model TEC500). 
The laser frequency is actively locked to the 4-3 hyper- 
fine component of the Cs Di transition (A = 894 nm) in 
an auxiliary cesium vapor cell by means of the dichroic 
atomic vapor laser lock (DAVLL) technique [IJ. The 
stabilization to a Doppler-broadened resonance provides 
a continuous stable operation over several weeks and 
makes the setup rather insensitive to mechanical shocks. 

At the sensor head the light from the fiber is coUi- 
mated by a / = 15 mm lens and its polarization is made 
circular by a polarizing beamsplitter and a quarter-wave 
plate placed before the cesium cell. The light transmitted 
through the cell is focused into the return fiber, which 
guides it to a photodiode. The photocurrent is ampli- 



fied by a low-noise transimpedance amplifier. Placing 
the laser, the electronics, and the photodiode far away 
from the sensor head eliminates magnetic interference 
generated by those components on the magnetometer (a 
photocurrent of 10 /iA, e.g., produces a magnetic field of 
200 pT at a distance of 1 cm) . In the present setup the 
oscillating-field coil is fed via a twisted-pair conductor, 
which represents an effective antenna by which electro- 
magnetic signals can be coupled into the magnetic shield. 
In a future stage of development it is planned to replace 
this electric lead by an opto-coupled system. 

Multimode fibers were used for ease of light coupling. 
We found that a few loops of 3 cm radius of curvature 
in the fiber led to quasi-depolarization of the initially 
linearly polarized beam, thereby suppressing noise con- 
tributions from polarization fluctuations. A rigid fixation 
of the flbers was found necessary to reduce power fluctu- 
ations of the flber transmission to a level of 4 x 10~^ in 
1 Hz bandwidth. 

The studies reported below were performed inside 
closed cylindrical shields consisting of three layers of 
MUMETALL (size of the innermost shield: length 
600 mm, diameter 300 mm) that reduces the influence 
of ambient magnetic field variations. For the measure- 
ment of the noise spectrum (Sec. IIVA|I and the study of 
the magnetic field stability (Sec. IIVG|I the shield was 
improved by three additional cylinders of CO-NETIC 
mounted inside of the MUMETALL shield (innermost 
diameter 230 mm). The longitudinal bias field of 2/iT, 
corresponding to a Cs Larmor frequency of 7 kHz, is pro- 
duced by a solenoid (length 600 mm, diameter 110 mm) 
inside the shield and the 8 mA current is provided by a 
specially designed stable current supply. 



A. Resonance linewidth 

The lineshapes of the magnetic resonance line are mea- 
sured with the magnetometer operating in the open-loop 
mode (Fig.Q] mode B). A sinusoidally oscillating current 
of frequency uj^f is supplied to the r.f. coils by a function 
generator, whose frequency is ramped across the Larmor 
frequency, and the output of the photodiode is demodu- 
lated by a lock-in amplifier. Magnetic resonance lines 
were recorded for different Bi amplitudes and different 
values of the pump light power. Typical resonance lines 
are shown in Fig. [3 The lineshapes were fitted by the 
function Q to the experimental (/'(I'rf) curves, which al- 
lows one to infer the linewidth AfnwHM- We recall that 
the linewidth is not affected by r.f. power broadening, 
but that it is subject to broadening by the optical pump- 
ing process. The dependence of Az^hwhm on the laser 
intensity (Fig. ^ shows that the optical broadening has 
a nonlinear dependence on the light intensity. The mini- 
mum or intrinsic linewidth is determined by extrapolat- 
ing A^-HWHM to zero light intensity. 

For a J = 1/2 two-level system theory predicts a lin- 
ear dependence of the linewidth on the pumping light 
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FIG. 2: Magnetic resonance spectra obtained by scanning 
the frequency Vri of the oscillating field: a) quadrature com- 
ponent, b) in-phase component, c) phase between the oscillat- 
ing field and the modulation of the transmitted power. The 
Larmor frequency vl is 7002.3 Hz, the power-broadened half 
linewidth is 2.2 Hz. The intrinsic half linewidth of 1.4 Hz is 
indicated. 



intensity, as long as stimulated emission processes from 
the excited state can be neglected. However, the mag- 
netic resonance spectrum in the F = 4 manifold of the 
Cs ground state is a superposition of eight degenerate 
resonances corresponding to all allowed AM = ±1 tran- 
sitions between adjacent Zeeman levels. The coupling of 
the CT+ polarized light to the different sublevels depends 
on their magnetic quantum number Mp and is given by 
the corresponding electric dipole transition matrix ele- 
ments. As a consequence each of the eight resonances 
broadens at a different rate. The observed linewidth 
results from the superposition of those individual lines 
weighted by the population differences of the levels cou- 
pled by the r.f. transition and the corresponding mag- 
netic dipole transition rates. The observed nonlinear de- 
pendence of the width on the light intensity follows from 
the nonlinear way in which those population differences 
and hence the relative weights are changed by the optical 
pumping process. 



We have calculated the lineshapes of the magnetic res- 
onance lines by numerically solving the Liouville equa- 
tion for the ground state density matrix. Interactions 
with the optical field as well as the static and oscillating 
magnetic fields were taken into account in the rotating 
wave approximation. We further assumed an isotropic 
relaxation of the spin coherence at a rate given by the 
experimentally determined intrinsic linewidth of Fig. |31 
The solid curve in that figure represents the linewidths 
inferred from the calculated lineshapes. The calculations 
used as a variable an optical pumping rate (proportional 
to the light power intensity) and the only parameter used 
to fit the calculation to the experimental data was the 
proportionality constant between the laser intensity and 
that pump rate. 

The intrinsic linewidth, i.e., the linewidth for vanish- 
ing optical and r.f. power, is determined by relaxation 
due to spin exchange Cs-Cs collisions, Cs-wall collisions, 
and collisions of the atoms with the Cs droplet in the 
reservoir sidearm. The latter contribution depends on 
the ratio of the cross section of the constriction in the 
sidearm and the inner surface of the spherical cell. With 
an inner sidearm diameter of 0.5 mm that contribution to 
the HWHM linewidth can be estimated to be on the or- 
der of Az^ — IHz. The contribution from spin exchange 
processes at room temperature to the linewidth can be 
estimated using the cross-section reported in [13| to be 
on the order of 3 Hz, which is larger than the measured 
width. A possible explanation for this discrepancy is the 
adsorption of Cs atoms in the paraffin coating 14], which 
may lead to an effective vapor pressure in the cell below 
its thermal equilibrium value. 



B. Magnetometer mode 

The actual magnetometry is performed in the phase- 
stabilized mode (Fig. 2] mode A) as described above. The 
photodiodc signal is demodulated by a lock-in amplifier 
(Stanford Research Systems, model SR830) locked to the 
driving r.f. frequency, produced by a voltage-controlled 
oscillator (VCO). The time constant of the lock-in am- 
plifier was set to r = 30 /iS, which corresponds to a band- 
width of 2.6 kHz with a -24dB/octave filter roll-off. Ei- 
ther the phase (adjusted to be 0° on resonance) or the 
dispersive in-phase signal of the lock-in amplifier can be 
used to control the VCO, and hence to lock its frequency 
to the center of the magnetic resonance. Compared to the 
in-phase signal the phase signal of the lock-in amplifier 
has the advantage that the resonance linewidth is not af- 
fected by r.f. power broadening. However, the bandwidth 
of the phase output of the digital lock-in amplifier used 
was limited to 200 Hz by its relatively slow update rate. 
For the neutron EDM experiment the magnetometer has 
to be operated with the highest possible bandwidth. We 
therefore chose the in-phase signal for the following stud- 
ies. That signal drives the VCO via a feedback amplifier 
(integrating and differentiating), which closes the feed- 
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Light intensity (|^W/mm^) 

FIG. 3: Resonance HWHM linewidth as a function of the 
light intensity delivered to the cell. The power of the laser 
beam is given by II ■ 2.8 mm^. The dots represent the widths 
obtained from the phase signal of the lock-in amplifier with 
very low r.f. power. The extrapolated intrinsic linewidth is 
1.4(1) Hz. The solid line is a one-parameter fit of a numerical 
calculation to the data (see text). The size of the symbols 
represents the vertical error bars. 



back loop locking the radio frequency to the Larmor fre- 
quency. 
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FIG. 4: Square root of the power spectral density (PSD) of 
the magnetometer output frequency relative to the Larmor 
frequency of uo — 6998 Hz (averaged 20 times). The lowest 
lying curve is the contribution from laser power noise and A'lnt 
its contribution at the Larmor frequency (inset). The straight 
line indicates the shot noise limit Nsn- The signal-to-noise 
ratio S/Nint is approximately 33000. The sidebands are due 
to imperfectly shielded magnetic field components oscillating 
at the 50-Hz power line frequency. The signal-to-noise ra- 
tio S'/A'ext due to external field fiuctuations is approximately 
4600. All measurements were performed in a 1 Hz bandwidth. 



IV. PERFORMANCE OF THE 
MAGNETOMETER 

A. Magnetometric sensitivity 

We characterize the sensitivity of the magnetometer 
in terms of the noise equivalent magnetic flux density 
(NEM), which is the flux density change 6B equivalent 
to the total noise of the detector signal 



SB^ ^ 5Bl 



5Bl 



(4) 



with both internal and external contributions: (5i?int de- 
scribes limitations due to noise sources inherent to the 
magnetometer proper, while (5Z?oxt represents magnetic 
noise due to external field fluctuations. The larger of 
the two contributions determines the smallest flux den- 
sity change that the magnetometer can detect. In general 
the internal NEM ffiint may have several contributions, 
which may be expressed as 



5BL = 
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Al^HWHM \ ' 
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At^HWHM \ 



(5) 
where S is the magnetometer signal, Nq^-^ are the noise 
levels of the different processes contributing to SBi^t , and 
iVgN the fundamental shot noise limit of the 0PM signal. 
7 is approximately 3.5Hz/nT for ^^■^Cs and Ai^hwhm is 
the half width of the resonance (cf. Sec lIVB)) . 
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The external NEM (5i?oxt can be parametrized in the 
form of Eq. (0) with an equivalent signal noise iVcxt so 
that Eq. I0J can be expressed as 



(6) 



In a strict sense Eqs. ^ and © are valid for the open 
loop operation of the magnetometer. The parameters 7 
and Ai^HWHM do not depend on the mode of operation, 
whereas S/N may very well be affected by feedback. If 
the bandwidth, in which the noise is measured, is much 
smaller than the loop bandwidth of 1 kHz (cf. Sec. lIVB|l . 
it is reasonable to assume that the power spectral density, 
i.e. S/N, is not affected by the loop being opened or 
closed. 

Experimentally the noise levels Na are determined 
from a Fourier analysis of the photodiode signal, when 
the magnetometer is operated in the phase-stabilized 
mode under optimized parameter conditions. Each noise 
level N is defined as the square root of the integrated 
(frequency dependent) power spectral density p| of the 
corresponding signal fluctuations 



1/2 



N 



Pl'^f 



(7) 



where /bw is the measurement bandwidth. If the noise 
is white or if the bandwidth is much smaller than the 



width of typical spectral features in the power spectrum 
the noise level at a given frequency / is given by 



N = psVfhw = ps/V2t, 



(8) 



where r is the integration time used for calculating the 
Allan standard deviation introduced below. Figure 01 
shows a typical Fourier spectrum of the 0PM signal. 
The prominent central feature corresponds to the Lar- 
mor oscillation of the photocurrcnt at 7 kHz during the 
phase-stabilized operation of the 0PM. It is the signal-to- 
noise ratio at the Larmor frequency which determines the 
NEM of the magnetometer. The Larmor peak (carrier) 
is superposed on a 20 Hz broad pedestal, which itself lies 
above a constant noise floor. The two discrete sidebands 
originate from magnetic fields oscillating at the 50 Hz 
line frequency. The pedestal results mainly from imper- 
fectly shielded low- frequency field fluctuations. The con- 
tinuous spectrum of such fluctuations is mixed with the 
Larmor frequency and therefore appears as a symmetric 
background underlying the carrier. We have fitted the 
pedestal with a Lorentzian lineshape, from which we in- 
fer a signal-to-noise ratio S/Ncxt = 4600 at the Larmor 
frequency, in which A'oxt corresponds to magnetic field 
fluctuations (5Bcxt of 210 fT in a 1 Hz bandwidth. 

A major contribution to the internal noise A'int comes 
from laser power fluctuations. This can be understood 
as follows. The detected photocurrcnt has the general 
form /pD — Ipc + A/cos(a'rf -I- 0), where /pc is the DC 
component and A/ the amplitude of the AC component, 
which depends on the atomic number density and the 
degree of spin polarization. If the magnetic resonance 
signal A/ contains fluctuations at frequency w, these will 
appear in /po at the frequencies jwrf ±a;|. Low frequency 
power will therefore produce a symmetric background 
under the Larmor peak, which also contributes to the 
pedestal discussed above. In order to distinguish it from 
the SBext contributions we have measured the low fre- 
quency noise spectrum of the laser power. It shows a 
l//-like behavior for small frequencies, which levels off 
at the shot noise value. After mixing with the Larmor 
frequency its contribution to the magnetometer noise is 
shown as the lowest lying curve in Fig. ^ The central 
portion of this spectrum is shown on an expanded scale 
in the inset. The sidebands are probably due to mechani- 
cal vibrations, while the central feature extends ±1.5 Hz. 
-^int represents the contribution of that power noise to 
the magnetometer noise at the Larmor frequency. The 
direct contribution to SB from light power noise is thus 
negligible compared to the contributions from field fiuc- 
tuations and confirms the assumption made above that 
the latter dominate the pedestal. The signal-to-noise ra- 
tio 5'/iVint is 33000 and yields a NEM SBint = 29 fT in a 
1 Hz bandwidth, determined by laser power fluctuations. 

Far away from the Larmor frequency the measured 
constant noise floor exceeds the calculated shot noise 
level (iVsN) by a factor of 1.5. This may originate from 
additional noise sources related, e.g., to the laser fre- 
quency stabilization. The fundamental limit of the mag- 



netometric sensitivity is determined by the white shot 
noise Nsn — \/2eIpcfhw of the DC component of the 
photocurrent, /pc, which deflnes the ultimate shot noise 
limited NEM (5Bsn- Under optimized conditions the pho- 
tocurrcnt is 5 /iA, which yields a shot noise limited NEM 
of SBsN = 10 fT in a bandwidth of 1 Hz. 

Light shift fluctuations are an additional source of 
noise. Any fluctuations of the parameters causing a light 
shift (laser power and/or laser frequency detuning) will 
produce magnetic fleld equivalent noise. We will show 
later that for a 1 Hz detection bandwidth this effect gives 
a negligible contribution to the Fourier spectrum. 

As the internal noise level SBi^t is much smaller than 
the external fleld fluctuations SB^xt the magnetometer is 
well suited to measure the characteristics of such fleld 
fluctuations (cf. Sec. IIVC|) and/or to compensate them 
using an active feedback loop. The accuracy of such mea- 
surements or the performance of such a stabilization is 
ultimately limited by the internal noise of the magne- 
tometer, which under ideal conditions can reach the shot 
noise limit. 



B. Magnetometer optimization and response 
bandwidth 



According to Eq. ^ the sensitivity of the magnetome- 
ter depends on the resonance linewidth Az^hwhm and on 
the signal-to-noise ratio. For given properties of the sen- 
sor medium (cesium vapor pressure and cell size) these 
two properties depend on the two main system param- 
eters, viz., the laser intensity II (or power Pl) and the 
amplitude Bi of the r.f. fleld. For the application in 
the neutron EDM experiment the sensor size and va- 
por pressure are dictated by the experimental constraints 
(flxed geometry and operation at room temperature), so 
that the experimental optimization of the magnetomet- 
ric sensitivity is performed in the (7^, Bi) space by an 
iterative procedure. Fig. |S1 shows examples of signal- 
to-noise ratio recordings during such an iteration. The 
optimum operating point was found for a laser inten- 
sity II of 9^W/mm and a r.f. field amplitude Bi of 
2.7 nT. The resonance linewidth under optimum condi- 
tions is Ai/HWHM = 3.4(1) Hz, which exceeds the intrinsic 
linewidth by a factor of 2.4. 

In order to investigate the dependence of the NEM 
on the two optimization parameters we have calculated 
that dependence using the density matrix formalism by 
assuming that the signal noise is determined by the shot 
noise of the photocurrent. The result is shown in Fig.|H|as 
a density plot. One recognizes a broad global minimum 
which is rather insensitive to the parameter values as it 
rises only by 5% when the optimum light and r.f. power 
are varied by 50%. 

The bandwidth of the magnetometer, i.e., its temporal 
response to field changes was measured in the following 
way: a sinusoidal modulation of the static magnetic field 
with an amplitude of 5 nT was applied by an additional 
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FIG. 5: Dependence of the experimental signal-to-noise ratio 
(measured in a 1 Hz bandwidth) on the light intensity and the 
current applied to the r.f. coils. The noise was measured 70 Hz 
away from the carrier. The r.f. current in a) was S^App, the 
light intensity in b) was 7/xW/mm'^. The dashed lines are 
drawn to guide the eyes. These are typical recordings used to 
optimize the system parameters. 



single wire loop (110 mm diameter) wound around the Cs 
cell. The response of the magnetometer to that perturba- 
tion was measured directly on the VCO input voltage in 
the phase-stabilized mode. The result is shown in Fig.[7| 
The overall magnetometer response follows the behavior 
of a low-pass filter (— 24dB/octave roll-off) with a -3dB 
point at approximately IkHz. The lock-in time constant 
was 30 /iS which corresponds to a bandwidth of 2.6 kHz. 
The difference is due to additional filters in the feedback 
loop. 



Application: Field fluctuations in a magnetic 
shield 



External field fluctuations are the dominant contribu- 
tion to the noise of the LsOPM when it is operated in 
the six-layer magnetic shield. We have used the mag- 
netometer to study the temporal characteristics of the 
residual field variations. The Allan standard deviation 




II (arb. units) 

FIG. 6: Density plot (in arbitrary units) of the calculated 
NEM SBsN, in dependence on the laser intensity II and the 
r.f. amplitude Bi. 
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FIG. 7: Frequency dependence of the magnetometer response 
to a small amplitude sinusoidal modulation of the static field 
-Bo (circles). The solid line indicates the amplitude response 
of a 4"^-order low-pass filter (-24 dB/octave roll-off). 



[l5| is the most convenient measure for that characteri- 
zation. With respect to the experimental specifications 
of the neutron EDM experiment our particular interest 
is the field stability for integration times in the range of 
100 to 1000 s. For that purpose we recorded the Larmor 
frequency in multiple time series of several hours with 
a sampling rate of 0.1s by feeding the photodiodc sig- 
nal, filtered by a resonant amplifier (band-pass of 200 Hz 
width centered at 7 kHz), to a frequency counter (Stan- 
ford Research Systems, model SR620). From each time 
series the Allan standard deviation of the flux density in- 
side the shield was calculated. A typical result is shown 
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FIG. 8: a) Allan standard deviation (SBcxt of the magnetic flux 
density inside the magnetic shield (•). b) NEM (JBint (□) lim- 
ited by laser power fluctuations; extrapolated NEM SBsn (B) 
for shot noise limited operation. The slope represents Eq. © 
assuming a white noise behavior, c) Measured contributions 
to SBint from light power fluctuations (t) with present setup. 
Solid lines in a) and c) are drawn to guide the eye. The dwell 
time of the frequency counter was 100 ms. 



in Fig.|Slwith both absolute and relative scales. For inte- 
gration times below one second the observed fluctuations 
(curve a) decrease as t~^/^, indicating the presence of 
white field-amplitude noise. It is characterized by a spec- 
tral density of 245 fT/VHz. Although the Allan standard 
deviation represents a different property than the Fourier 
noise spectrum it is worthwhile to note that the latter 
value is comparable with the NEM (Si?oxt = 210 fP of the 
pedestal in Fig. 0] discussed above. The field fiuctuations 
reach a minimum value of approximately 240 fT for an 
integration time of 0.7 s. 

The central region of the Allan plot (Fig. |Ht) shows a 
bump for integration times of 1-100 seconds. It is prob- 
ably due to fiuctuations of the 8 mA current producing 
the 2//T bias field. A magnetic field fiuctuation of 200 fT 
corresponds to a relative current stability of lO^'', i.e., 
to current fiuctuations of 800 pA. In an auxiliary experi- 
ment we measured the current fiuctuations A/ by record- 
ing voltage fiuctuations over a series resistor for several 
hours. We found relative fluctuations of A/// in the cor- 
responding Allan plot of the same order of magnitude as 
the AB/B fluctuations. It is thus reasonable to assume 
that the origin of the plateau in Fig.|SjL is due to current 
fiuctuations of the power supply. 

The measurement of the magnetic field during several 
days shows fiuctuations with a period of one day and 
an amplitude of about 1 Hz, superposed by additional 
uncorrelated drifts. The periodic fiuctuations are proba- 
bly due to changes of the solenoid geometry induced by 
temperature fiuctuations. The Allan standard deviations 
for integration times exceeding 200 s are thus determined 
by temperature fiuctuations and drifts of the laboratory 
fields, which are not completely suppressed by the shield. 



D. Frequency noise due to light power fiuctuations 

It is well-known that a near-resonant circularly polar- 
ized light field shifts the Zeeman levels in the same way 
as a static magnetic field oriented along the light beam. 
The light shift has contributions from the AC Stark shift 
and coherence shift due to virtual and real transitions 
[ig |. The AC Stark shift, and hence the equivalent mag- 
netic field Bls is proportional to the light intensity II 
and has a dispersive (Lorentzian) dependence on the de- 
tuning of the laser frequency from the center of the op- 
tical absorption line. It is therefore expected to vanish 
at the (optical) line center. In our experiment the laser 
frequency is locked to the center of a Doppler-broadened 
hyperfine component. However, that frequency does not 
coincide with the frequency for which the light shift van- 
ishes, because of finite light shift contributions from the 
adjacent hyperfine component. While the two hyperfine 
components are well separated in the optical absorption 
spectra, their corresponding light shift spectra overlap 
because of the broad wings of their dispersive lineshapes. 

In order to measure the light shift effect we period- 
ically changed the light power between P + AP/2 and 
P — AP/2 and recorded the corresponding Larmor fre- 
quencies. False effects from drifts of the external mag- 
netic field were suppressed by recording data over several 
modulation periods. For each modulation amplitude AP 
the Larmor frequency was measured with both (T+ and 
(T~ polarizations by rotating the quarter-wave plate by 
means of a mechanical remote control from outside the 
shield. 

The induced changes of the magnetometer readings for 
both polarizations are shown in Fig. ^ As anticipated, 
the shift of the Larmor frequency is proportional to the 
modulation amplitude of the light power and changes sign 
upon reversing the light helicity. However, it can be seen 
that the slope of the light shift depends on the helicity. 
This asymmetry is the result of contributions from three 
distinct effects, which we discuss only qualitatively here. 

(1) The light shift due to virtual transitions (AC Stark 
shift), which is proportional to the helicity of the light 
and thus leads to a symmetric contribution to the curves 
of Fig. 1^1 (equal in magnitude, but opposite in sign); (2) 
the light shift due to real transitions (coherence shift) [I^, 
whose origin is a change of the effective g-factor of the 
Cs atom due to the fact that with increasing laser power 
the atom spends an increasing fraction of its time in the 
excited state with a 3 times smaller gj;--factor of opposite 
sign than that of the ground state; (3) a possible power 
dependent change of the capacity of the photodiode and 
a subsequent power dependent phase shift of the pho- 
tocurrent. The latter two effects yield shifts which have 
the same sign for both light polarizations, so that the 
combined contribution of the three effects may explain 
the different magnitudes of the slopes. A quantitative 
study of those effects is underway. 

Using curve (a) as a worst-case estimate for the fiuctu- 
ations of the Larmor frequency due to light power fiuc- 
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netometer has an intrinsic NEM of 29 iT, defined as the 
Allan standard deviation for an bandwidth of IHz. If 
the 1// noise of the laser power can be lowered, e.g., by 
an active power stabilization and the excess white noise 
floor can be reduced to the shot noise level the LsOPM 
should reach a NEM of 10 fT for a 1 Hz bandwidth. The 
bandwidth of the phase-stabilized LsOPM is 1 kHz. We 
have used the LsOPM to measure field fluctuations in a 
six-layer magnetic shield for integration times between 
0.1 and 1000 seconds, whose lowest values were found to 
be on the order of 200-300 fT. Light shift fluctuations, 
against which no particular precautions were taken, are 
one to two orders of magnitude smaller than the residual 
field fluctuations in the shield. 



FIG. 9: Relative light shift of the Larmor frequency as a 
function of the relative modulation amplitude AP of the laser 
power P. Curves (a) and (b) represent measurements with 
light of opposite circular polarization. 



tuations we estimated, based on measured power fluc- 
tuations, the resulting magnetic field fiuctuations. The 
results are shown as triangles in Fig. |S1 Light shift fluc- 
tuations of the magnetometer readings are thus one to 
two orders of magnitude smaller than residual field fluc- 
tuations in the present shield. The light shift noise can 
of course be further suppressed by adjusting the laser 
frequency to the zero light shift frequency point or bet- 
ter by actively stabilizing it to that point or by actively 
stabilizing the laser power. 



V. SUMMARY AND CONCLUSION 

We have described the design and performance of a 
phase-stabilized cesium vapor magnetometer. The mag- 



The LsOPM described here compares very favorably 
with state-of-the-art lamp-pumped magnetometers. De- 
tails on that comparison will be published elsewhere. It 
will be a valuable tool for fundamental physics experi- 
ments. The LsOPM presented above meets the require- 
ments of the neutron-EDM experiment on the relevant 
time scales in the range of 100 to 1000 seconds. 
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